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ABSTRACT 

Recently, it was reported that an electrostrictive graft elastomer exhibits large electric field-induced strain (4%). Combined 
with its high mechanical modulus, the elastomer can offer very promising electromechanical properties, in terms of output 
mechanical energy density, for an electroactive polymeric material. Therefore, it has been considered as one of the 
candidates that can be used in high performance, low mass actuation devices in many aerospace applications. Various bi- 
layer-based bending actuators have been designed and fabricated. An analytic model based on beam theory in the strength of 
materials has been derived for the transverse deflection, or curvature, and the longitudinal strain of the bi-layer beam. The 
curvature and strain are functions of the applied voltage and the thickness, width, and Young's modulus of the active and 
passive layers. The model can be used to optimize the performance of electrostrictive graft elastomer-based actuators to meet 
the requirements of various applications. In this presentation, optimization and sensitivity studies are applied to the bending 
performance of such actuators. 
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1. INTRODUCTION 

Materials that sustain mechanical displacement under controlled electrical excitation are needed as actuators for many 
applications. For aerospace applications, there is also a need for low mass, high performance, and ease of processability, 
which are inherent characteristics of electroactive materials. Electroactive polymeric elastomers that show electromechanical 
activity, especially large electric field induced strain, are being considered as potential candidates for such applications. 
Existing materials include polyurethane elastomers [1-3] and silicone rubber [4,5]. Recently, we have demonstrated an 
electrostrictive response in graft elastomers [6,7]. The elastomers offer large electric field induced strain and high 
mechanical modulus. Therefore, a high electromechanical output power, or high strain energy density is achieved. In 
addition to the high performance as a new class of electromechanically active polymeric materials, the electrostrictive graft 
elastomers also offer such advantages as excellent processability and electrical and mechanical toughness. 

Smart materials, such as sensor and actuator materials, are finding many applications in space technologies. However, 
lightweight, low cost, and high performance are still challenges in the development of these types of functional materials. 
Recently, an electrostrictive graft elastomeric polymer was discovered at NASA Langley Research Center. It offers large 
electric field induced strain (4%). This material can be used in many NASA space technologies that include position and 
tension control of functional membranes, aerodynamic and flow control of spacecraft, and development of self-adjustable 
multifunctional membranes as well as micro-electro-mechanical-systems (MEMS). Using the material, several prototype 
actuators have been developed. It was found that the performance of the actuators using the electrostrictive graft elastomers 
depends on the configuration, materials selection, and fabrication. This work is intended to help optimize the performance of 
the actuators by means of design and modeling. 

2. BASIC CONFIGURATION 

Figure 1 shows the basic configuration of a bending actuator that uses electrostrictive graft elastomers. The active layer was 
coated with gold electrodes on opposite sides of the graft elastomer film. The active layer was bonded to a passive graft 
elastomer layer that has no electrode coating. The adhesive was a thin epoxy layer. Figures 2a and 2b show the fabricated 
electrostrictive graft elastomer-based bending actuator in unexcited and excited states, respectively. 
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Figure 1: Basic configuration of a bending actuator using the electrostrictive graft elastomers. 
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Figure 2: The bending actuator in the unexcited state (a) and the excited state (b). 


3. DERIVATION AND OPTIMIZATION OF MODEL 


3.1 Derivation 

We shall attempt to develop an analytic model for the transverse curvature K and longitudinal strain £ of a two- layer, 
unrestrained beam with one layer electrostrictive and the other passive. The model is based on beam theory in the strength of 
materials, as presented in [8]. Of the two layers comprising the beam, the electrostrictive layer is layer 1 and the passive 
layer is layer 2. The unexcited model beam is shown in Figure 3. In the unexcited state, the beam is flat and of length L . 
The height h , width b , and Young's modulus E can be different for the two layers. Plane cross-sections of the beam 
are assumed to remain plane under flexure so that strains vary directly with their distance from the neutral axis, where the 
flexural strain is zero. Strains are assumed not to exceed the yield point. 
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Figure 3: Model of a two-layer flat beam of length L , 
width b , and layer thickness h . The longitudinal 
direction lies along the y - axis, and the transverse direction 
along the Z - axis. 


For the laboratory prototype shown in Figure 2, L = 0.035 m, h { = h 2 =21x10 6 m, and /? = b, =5.75x10 3 m. 
The appropriate equations for developing this model are presented by F. L. Singer [8]. As in the treatment of thermal 
stresses, the electrostrictive layer is allowed to freely strain under excitation. To maintain longitudinal force equilibrium, the 
electrostrictive layer is then brought into compression and the passive layer into tension. The resulting longitudinal strain £ 
of the beam is given by 


£ = 


f E 2 b 2 h 2 | ^ 

v E AK , 


(la) 


where e a is the free strain of the electrostrictive layer, as given by 


<? = ml 


' ] O 2 

yK j 


(lb) 


where m = 1.535x10 18 m 2 (volt)" 2 is the electrostrictive constant for the graft elastomer, and V l is the voltage applied 
across the electrostrictive layer. With ( lb) substituted for e a , the longitudinal strain of the beam becomes 


£ = 


mV { 


f E 2 b 2 h 2 + 1 A 




J 


(lc) 
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The corresponding bending moment M is given by 


M = 


(A +h 2 )e a 


1 1 

- + - 


EA^h E 2 b 2 h 2 j 


( 2 ) 


The beam curvature K = 1/p , where p is the radius of curvature that results from this moment, is given by 


K = 


M 


( 3 ) 


where / is the moment of inertia about the neutral axis of the equivalent homogeneous beam having Young’s modulus 
E 1 everywhere. In this equivalent homogeneous beam, width /? is replaced by its equivalent width (/? E f )/E n . The 
neutral axis is the centroidal axis of the equivalent cross-section. The resultant expression for 1 is given by 


j _ {Efafh* +2E 1 b 1 E 2 b 2 h 1 h 2 feh 1 2 +3h 1 h 2 +2h 2 2 )+{E 2 b 2 ) 2 h 2 * 

1 2E 2 (E l b l h l +E 2 b 2 h 2 ) 

and the ciuvature K becomes 

K = 6E l b l E 2 b 2 h l h 2 (h ] +h 2 )e a 

{E l b l YhS + 2E x b x E^b^h x h^ ( 2 h { 2 + 3h l h 2 + 2 hp )+ [E^b n )~h 

Substituting ( lb) for e a and dividing the numerator and denominator by E ] b ] E J) n give finally 


(4) 


(5a) 


K = 


6 h 2 (h l +h 2 )mV l 


h x ah* + 2h [ h 2 ( 2 hp + 3h { h 2 + 2h 2 )+ 


a 


where 


a = 


EA 

E 2 b 2 


(5b) 


(5c) 


This expression is a general formula for the unrestrained curvature of a two-layer beam with layer 1 electrostrictive and layer 
2 passive. 

3.2 Optimizations 

Besides m and V ] , the curvature K is a function of , h n , and (X . We can maximize the curvature with respect to 
a by taking c)k/c)(X = 0 , which yields the (X - optimization condition 


a = 


( 6 ) 
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When this condition is met, d K/d(X~ is negative, which confirms that K is a maximum with respect to (X . Therefore, 
for given values of h { and h 2 , condition (6) gives the value of (X that maximizes the curvature K . Also, when (6) is 
satisfied, the curvature K is much less sensitive to small changes in (X (small changes in E { , b ] , £ , , or b n ) than to 
small changes in h t or h 2 . 

We can alternatively maximize the curvature K with respect to h 2 by taking c) K/c)h n =0, which yields the h 2 - 
optimization condition 

2 h 2 +3h 2 h ] -a/ij 3 = 0 (7) 

This cubic equation has only one root for h n that is real and positive for practical values of (X 

h 2 = f(a)\ (8a) 


where 


/(«) = 


1 

2 


• — 1 + 2 cc + 2 (cc 




(8b) 


A plot of f{cc) versus (X 


is shown in Figure 4. 



Figure 4: Plot of f (cc) versus CX , where 
f(pc ) is given by (8b) and CC by (5c). 


Thus, for any given values of h x and CX , the value of h n that maximizes the curvature K is given by (8a). 
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If we attempt to maximize the curvature K with respect to h } by taking c)K / 3/ij = 0 , we find that there arc no positive 
real solutions for h } . Therefore, for given values of h 2 and CC , there is no value for h [ that maximizes K . The reason 

that no optimum value for h { exists is clarified by plotting contours of In (?c) on the h { — h 2 plane, as shown in 
Figure 5. 


xl0“ 6 



Figure 5: Contour plot of In ( K ) 
on the /tj — h , plane. Relevant 
parameters are listed in Table 1. 


xl0“ 6 


As h 2 increases along a vertical line, K increases until h 2 reaches the value given by (8a), after which K decreases. In 
contrast, as decreases along a horizontal line, K continues to increase monotonically without leveling off to a 

maximum value. Therefore, for given values of h x and CC , there exists a value of h n that maximizes K ; but for given 
values of h , and OC , one can always increase K further by decreasing . 

4 . SENSITIVITY STUDIES 


At this point the model is not sufficiently close to the laboratory prototype to make a quantitative comparison. The 
laboratory prototype actually contains five layers, as shown in Figure 1. Besides the electroactive graft elastomer layer and 
the passive graft elastomer layer of the model, there is the adhesive layer between them and two metallic gold layers that coat 
the electroactive layer and across which the voltage is applied. The gold layers reduce the effective longitudinal strain of the 
electroactive layer and increase the beam’s resistance to bending. The adhesive layer also has similar effects. Eventually we 
hope to bring the model and prototype into agreement. At present, we will only attempt to gain qualitative insights from the 
model by performing sensitivity studies. 
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The sensitivity studies apply to a simplified version of the laboratory prototype where the effects on curvature of the two gold 
layers and the adhesive layer are neglected. The first study starts from a baseline state (indicated by the superscript (0)) that 

is OC - optimized by satisfying (6). The applied voltage V, 0 =215 volts, which gives a curvature to the beam of 

K° =5.75 m' 1 and radius of curvature p° =0.174 m. The other baseline parameters for sensitivity study 1 are given in 
Table 1. 


Table 1. OC -Optimized Baseline State For Sensitivity Study 1 


E,° = E 2 = 580x10 6 Pa 
b t ° = b 2 ° = 5.75xl(T 3 m 
h t ° =h 2 = 21xl0” 6 m 
7° = 3.55xl0“ 17 m 4 
m° = 1.535 xlO -18 m 2 (volts)' 2 

y,° =215 volts 

e a ° = 1.61xl0” 4 
M° = 1.18x10 7 Nm 
K° =5.75 m' 1 
p° =0.174 m 
e° = 8.04 xlO” 5 
L° = 0.035003 m 


Because the baseline longitudinal strain £° due to the baseline voltage V t ° , as given in Table 1, is so small, the increase in 
the beam length L° is negligible, and it remains so throughout both sensitivity studies. 

In each sensitivity study, various design parameters are varied one at a time, and we calculate the effects on the unrestrained 
curvature K (5b) and longitudinal strain £ (lc), as shown in Table 2. 


Table 2. Sensitivity Study 1 For OC -Optimized Baseline State 



/ 0 
e a !e a 

M / M° 

7/7° 

K/K° 

£/£° 

m/m° =1.1 

1.1 

1.1 

1 

1.1 

1.1 

Vi /Vi° =1.1 

1.21 

1.21 

1 

1.21 

1.21 

hjh* = .9 

1.23 

1.11 

.86 

1.30 

1.17 

h 2 /h 2 =.9 

1 

.9 

.86 

1.05 

1.05 

EJE ° =1.1 

1 

1.05 

1.05 

1 

1.05 

E 2 /E 2 °=. 9 

1 

.95 

1.05 

1 

1.05 

bjb" =1.1 

1 

1.05 

1.05 

1 

1.05 

b 2 /b 2 ° = .9 

1 

.95 

.95 

1 

1.05 



From this table, we see that a 10 percent increase in the electrostrictive constant m results in a 10 percent increase in the 
curvature K and a 10 percent increase in the longitudinal strain £ of the beam. The increase in K results from a 10 
percent increase in the bending moment M due to a 10 percent increase in the free strain e a of the active layer. This 
increase in e a also accounts for the increase in £ . In the formulas for K and £ , the voltage V l enters as a squared 

quantity. Therefore, a 10 percent increase in V l results in a 21 percent increase in the curvature K and a 21 percent 
increase in the longitudinal strain £ . 

A 10 percent decrease in the thickness of the active layer h x gives a 30 percent increase in the curvature K and a 17 
percent increase in £ . In this case, the increase in K results from an 1 1 percent increase in the bending moment M and 
a 14 percent decrease in the moment of inertia / . A 10 percent decrease in the thickness of the passive layer h 2 gives a 
5 percent increase in K and a 5 percent increase in £ . These changes result from a 10 percent decrease in M and a 14 
percent decrease in 1 . We infer that the curvature K of a beam could be tailored efficiently by varying the thickness of 
the active layer along its length. However, for the (X - optimization to remain valid for large thickness changes, both 
layers may need to be tailored. 

Ten percent changes in Young's moduli E ] and E 1 or in the widths /? and /?, can alter the longitudinal strain £ by 
5 percent but have a negligible effect on the curvature K . The insensitivity of K to these changes confirms that 
dxldoc = 0 when the (X - optimization condition (6) is satisfied. 

The second sensitivity study starts from a baseline state that is h 2 - optimized. According to Figure 4, when CC = 1 , which 

corresponds to the laboratory prototype, J(CX) = 1/2 and h 2 = h t / 2 . This value for h 2 deviates from the laboratory 
prototype and from sensitivity study 1. The baseline state for the second sensitivity study is given in Table 3, where values 
that differ from the first study are indicated by arrows ( — >) . In particular, the baseline curvature k" is substantially 
greater in Table 3 than in Table 1. 

Table 3. h n -Optimized Baseline State For Sensitivity Study 2 


E,° =E 2 ° =580x10 6 Pa 
b t ° =b 2 =5.75xl0' 3 m 
h" = 21xl0” 6 m 
— > /z 2 ° = 10.5 xlO” 6 m 
-4 7° =1.5xl(T 17 m 4 

m° = 1.535xl0” 18 m 2 (volts)- 2 
Vj° =215 volts 
e a ° =1.61xl0“ 4 
— > M° =5.92xl<r 8 Nm 

-XK° =6.81 m' 1 

-4 P° =0.147 m 
— > £° = 1.07 xlO” 4 
-4 L° = 0.035004 m 
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As before, the baseline values are varied one at a time by 10 percent, and the effects on curvature K , longitudinal strain £ , 
and other relevant quantities are given in Table 4. 


Table 4. Sensitivity Study 2 For h 2 -Optimized Baseline State 



1 0 
e e 

a a 

M / M° 

III 0 

K/K° 

£/£° 

m/m° =1.1 

l.i 

1.1 

1 

1.1 

1.1 

VJV* =1.1 

1.21 

1.21 

1 

1.21 

1.21 

hjh* =: 9 

1.23 

1.11 

.81 

1.37 

1.19 

h 2 /h 2 ° = .9 

1 

.9 

.9 

1 

1.03 

EJE* =. 9 

1 

.96 

.95 

1.02 

.96 

e 2 /e 2 °= 1.1 

1 

1.06 

.95 

1.02 

.97 

bjb? = .9 

1 

.96 

.95 

1.02 

.96 

b 2 /b 2 ° =1.1 

1 

1.06 

1.05 

1.02 

.97 


A 10 percent decrease in h { gives a 37 percent increase in K . Therefore, besides greater baseline curvature, the h 2 - 
optimized state gives greater curvature sensitivity than the OC - optimized state. This result could change for different 
voltages, etc., and both optimizations should be considered in the design process. As expected for the h 2 -optimization, the 
curvature K is insensitive to small changes in h 2 ; however, K now has gained some sensitivity to E ] , b x , E 2 , and 
b 2 , which was missing in the (X - optimization. 

5. CONCLUDING REMARKS 

Beam theory in the strength of materials was used to obtain an analytic model for the transverse curvature and longitudinal 
strain of a two-layer actuator consisting of one electrostrictive layer and one passive layer. The layers can have different 
thicknesses, widths and Young's moduli. The curvature formula was optimized in two ways: (1) by varying the thickness of 
the passive layer and (2) by varying the ratio of the product of width and Young's modulus for both layers. The latter 
optimization was satisfied by a simplified version of the laboratory prototype. The former optimization required halving the 
thickness of the prototype’s passive layer. Two sensitivity studies (one for each optimization) were performed with the 
model. Both studies showed that a 10 percent increase in the voltage across the electrostrictive layer gives a 21 percent 
increase in the curvature. A 10 percent decrease in the thickness of the electrostrictive layer gives a 30 percent increase in 
curvature for the optimization satisfied by the laboratory prototype and a 37 percent increase for the other optimization. The 
curvature is much less sensitive to small variations in the thickness of the passive layer. Therefore, the curvature of an 
actuator beam can be tailored efficiently along its length by varying the thickness of its active layer. However, for either 
optimization to remain valid for large thickness changes, both layers may need to be varied. 
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